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Abstract 

Data from e"'"e~ annihilation into hadrons collected by the JADE experiment at 
centre-of-mass energies between 14 GeV and 44 GeV were used to study moments of 
event shape distributions. The data were compared with Monte Carlo models and 
with predictions from QCD NLO order calculations. The strong coupling constant 
measured from the moments is 

as(Mzo) = 0.1286 ± 0.0007(stat.) ± 0.0011(exp.) ± 0.0022(had.) ± 0.0068(theo.) , 
q;s(Mzo) = 0.1286 ± 0.0072 (total error) 

consistent with the world average. However, systematic deficiencies in the QCD 
NLO order predictions are visible for some of the higher moments. 



This note describes preliminary JADE results 



^See [1] for the full list of authors 



1 Introduction 



The annihilation of an electron and a positron into hadrons allows precise tests of Quan- 
tum Chromodynamics (QCD). Commonly jet production rates or distributions of event 
shape observables have been studied. Predictions by perturbative QCD combined with 
hadronization corrections derived from models have been found to describe the data at 
low and high energies well, see e.g. [2-4]. 

In this analysis we used data collected by the JADE experiment in the years 1979 to 
1986 at the PETRA e+e" coUider at DESY at six centre-of-mass energies ^/s covering the 
range of 14-44 GeV. We measured the first five moments of event shape observables for 
the first time in this ^/s region in e+c^ annihilation and compared the data to predictions 
by Monte Carlo Models and by perturbative QCD. From the comparison of the data with 
the theory we extracted the strong coupling constant as- 

The outhne of the note is as follows. In section 2, we define the observables used in 
the analysis and describe the perturbative QCD predictions. In section 3 the analysis 
procedure is explained in detail. Section 4 contains the discussion of the systematic 
checks which were performed and the resulting systematic errors. We collect our results 
in section 5 and summarize them in section 6. 



The properties of hadronic events may be characterized by a set of event shape observables. 
These may be used to characterize the distribution of particles and thus the topology of 
an event as "pencil-like", planar, spherical etc. They can be computed either using the 
measured charged particle tracks and calorimeter clusters, or using the true hadrons or 
partous in simulated events. The following event shapes were considered here: 

Thrust T: defined by the expression [5,6] 



The thrust axis Ut is the direction n which maximizes the expression in parentheses. 
A plane through the origin and perpendicular to tit divides the event into two 
hemispheres Hi and H2. 

C-parameter: The linearized momentum tensor 6"^ is defined by 



Heavy Jet Mass M^: The hemisphere invariant masses are calculated using the par- 
ticles in the two hemispheres Hi and H2. We define Mh [8,9] as the heavier mass, 
divided by \/s. 



2 Observables 





The three eigenvalues \j of this tensor define C [7] through 



C = 3(AiA2 + A2A3 + A3A1) . 
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Jet Broadening observables and -Bw^ These are defined by computing the quan- 
tity 

for each of the two event hemispheres, H^, defined above. The two observables [10] 
are defined by 

Bt — Bi + B2 , and = max(i?i, B2) 
where B^ is the total and i?w is the wide jet broadening. 

Transition value between 2 and 3 jets 2/23: Jet algorithms arc applied to cluster the 
large number of particles of an hadronic event into a small number of jets, reflecting 
the parton structure of the event. For this analysis we used the Durham scheme [11]. 
Defining each particle initially to be a jet, a resolution variable is calculated for 
each pair of jets i and j: 

Vij = ^^2 - COS , (1) 



where Ei and Ej are the energies, cosOij is the angle between the two jets and 
£^vis is the sum of the energies of all visible particles in the event (or the partons 
in a theoretical calculation). If the smallest value of yij is less than a predefined 
value Ucut, the pair is replaced by a jet with four momentum pfj — pf +Pj, and 
the clustering starts again with instead of the momenta pf and pj. Clustering 
ends when the smallest value of yij is larger than i/cut- The remaining jets arc then 
counted. The value of y^ut at which for an event the transition between a 2-jet and 
a 3-jet assignment occurs is called 7/23- 

In the following discussion, whenever we wish to refer to a generic event shape observ- 
able we use the symbol y. In almost all cases, larger values of y indicate regions dominated 
by the radiation of hard gluons and small values of y indicate the region infiuenced by 
multiple soft gluon radiation. Note that thrust T forms an exception to this rule, as the 
value of T reaches unity for events consisting of two coUimated back-to-back jets. Wc 
therefore use y = 1 — T instead. For all of these event shapes, a perfectly coUimated 
("pencil-like") two-jet final state will have y = 0. 

The nth, n = 1,2,..., moment of the distribution of an event shape observable y is 
defined by 

(!/") = / y''-:r^y , 

JO a ay 

where ymax is the kinematically allowed upper limit of the observable. The calculations 
always involve a full integration over phase space, which implies that comparison with 
data always probes all of the available phase space. This is in contrast to QCD predictions 
for distributions; these are commonly only compared with data, e.g. in order to measure 
as- in restricted regions, where the theory is able to describe the data well, see e.g. [2]. 
Comparisons of QCD predictions for moments of event shape distributions with data are 
thus complementary to tests of the theory using the differential distributions. 
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The formula for the 0{al) QCD prediction of (y") is with as = Q;s/(27r) 

(l/") = Anas + Bncxl (2) 

involving the 0{as) coefficients An and 0(as) coefficients Bn- The values of the coeffi- 
cients An and Bn can be obtained by numerical integration of the QCD matrix elements 
using the program EVENT2 [12]. 

The QCD prediction depends on the renormahzation scale ^, see e.g. [13]. The renor- 
malization scale factor is defined as = H-j \fs implying that as = as(/^) in equ. (2). A 
truncated fixed order QCD calculation such as equ. (2) will depend on x^. The renormal- 
ization scale dependence is implemented by the replacement i3„ — > i3„ + /3o In x^An where 
/3o = 11 — \ni is the leading order /3- function coefficient of the renormahzation group 
equation and rif = 5 is the number of active quark flavours. 

The data are normalized to the total hadronic cross section crtot — cro(l + Q;s/7r) in LO 
while the QCD calculations are normalized to the born cross section ctq. Thus a correction 
has to be applied by making the replacement Bn ^ Bn — 2An- 

3 Analysis Procedure 

3.1 The JADE Detector 

A detailed description of the JADE detector can be found in [1]. This analysis relies mainly 
on the reconstruction of charged particle trajectories and on the measurement of energy 
deposited in the electromagnetic calorimeter. Tracking of charged particles was performed 
with the central detector, which was positioned in a solenoidal magnet providing an axial 
magnetic field of 0.48 T. The central detector contained a large volume jet chamber. Later 
a vertex chamber close to the interaction point and surrounding ^-chambers to measure 
the 2;-coordinate ^ were added. Most of the tracking information was obtained from the 
jet chamber, which provided up to 48 measured space points per track, and good tracking 
efficiency in the region | cos^| < 0.97. Electromagnetic energy was measured by the lead 
glass calorimeter surrounding the magnet coil, separated into a barrel (| cos < 0.839) and 
two end-cap (0.86 < |cos^^| < 0.97) sections. The electromagnetic calorimeter consisted 
of 2520 lead glass blocks with a depth of 12.5 radiation lengths in the barrel (since 1983 
increased to 15.7 in the middle 20% of the barrel) and 192 lead glass blocks with 9.6 
radiation lengths in the end-caps. 

3.2 Data Samples 

The data used in this analysis were collected by JADE between 1979 and 1986 and corre- 
spond to a total integrated luminosity of 195 pb^^. The breakdown of the data samples, 
mean centre-of-mass energy, energy range, data taking period, collected integrated lumi- 
nosities and the size of the data samples after selection of hadronic events are given in 

^In the JADE right-handed coordinate system the +x axis pointed towards the centre of the PETRA 
ring, the y axis pointed upwards and the z axis pointed in the direction of the electron beam. The polar 
angle 6 and the azimuthal angle (j) were defined with respect to z and x, respectively, while r was the 
distance from the z-axis. 
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table 1. The data samples were chosen following previous analyses, e.g. [1-4,14]. The 

data are available from two versions of the reconstruction software from 9/87 and from 
5/88. We used both sets and considered differences between the results as an experimental 
systematic uncertainty. 



average 
energy in GeV 


energy 
range in GeV 


year 


luminosity 
(Pb-^) 


selected 
events 

9/87 


selected 
events 

5/88 


14.0 


13.0-15.0 


1981 


1.46 


1722 


1783 


22.0 


21.0-23.0 


1981 


2.41 


1383 


1403 


34.6 
35.0 


33.8-36.0 
34.0-36.0 


1981-1982 
1986 


61.7 
92.3 


14213 
20647 


14313 
20876 


38.3 


37.3-39.3 


1985 


8.28 


1584 


1585 


43.8 


43.4-46.4 


1984-1985 


28.8 


3896 


4376 



Table 1: The average center-of-mass energy, energy range, year of data taking and inte- 
grated luminosity for each data sample, together with the numbers of selected data events 
using the data sample version of 9/87 or 5/88. 



3.3 Monte Carlo Samples 

Samples of Monte Carlo simulated events were used to correct the data for experimen- 
tal acceptance and backgrounds. The process e'^'e" — > hadrons was simulated using 
PYTHIA 5.7 [15]. Corresponding samples using HERWIG 5.9 [16,17] were used for sys- 
tematic checks. The Monte Carlo samples generated at each energy point were processed 
through a full simulation of the JADE detector [18-20], summarized in [14], and recon- 
structed in essentially the same way as the data. 

In addition, for comparisons with the corrected data, and when correcting for the 
effects of fragmentation, large samples of Monte Carlo events without detector simula- 
tion were employed, using the parton shower models PYTHIA 6.158, HERWIC 6.2 and 
ARIADNE 4.11 [21]. All models were adjusted to LEP 1 data by the OPAL collaboration. 

The ARIADNE Monte Carlo generator is based on a color dipole mechanism for the 
parton shower. The PYTHIA and HERWIG Monte Carlo programs use the leading 
logarithmic approximation (LLA) approach to model the emission of gluons in the parton 
shower. For the emission of the first hard gluon, the differences between the LLA approach 
and a leading order matrix element calculation are accounted for. However, the emission 
of gluons later on in the parton shower is based only on a LLA cascade and is expected to 
differ from a complete matrix element calculation. For this reason we do expect deviations 
in the description of the data by the Monte Carlo models. Recently new Monte Carlo 
generators have been developed, which implement a more complete simulation of the hard 
parton emission [22]. However, the models are not yet tuned to data taken at LEP and 
were therefore not considered in this analysis. 
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3.4 Selection of Events 



The selection of events for this analysis aims to identify hadronic event candidates and 
to reject events with a large amount of energy emitted by initial state radiation (ISR). 
The selection of hadronic events was based on cuts on event multiplicity (to remove 
leptonic final states) and on visible energy and longitudinal momentum balance (to remove 
radiative and two-photon events, e+e" —>■ e+e" hadrons). The cuts used are documented 
in [23-25] and summarized in a previous publication [2] . 

Standard criteria were used to select good tracks and clusters of energy deposits in 
the calorimeter for subsequent analysis. Charged particle tracks were required to have at 
least 20 hits in r-0 and at least 12 in r-z in the jet chamber. The total momentum was 
required to be at least 50 MeV. Furthermore, the point of closest approach of the track 
to the collision axis was required to be less than 5 cm from the nominal collision point in 
the X — y plane and less than 35 cm in the direction. 

In order to mitigate the effects of double counting of energy from tracks and calorime- 
ter clusters a standard algorithm was adopted which associated charged particles with 
calorimeter clusters, and subtracted the estimated contribution of the charged particles 
from the cluster energy. Charged particle tracks were assumed to be pions while the 
photon hypothesis was assigned to electromagnetic energy clusters. Clusters in the elec- 
tromagnetic calorimeter were required to have an energy exceeding 0.15 GeV after the sub- 
traction of the expected energy deposit of any associated tracks. From all accepted tracks 
and clusters i the visible energy E^is = J2iEi, momentum balance pbai = \ J2iPz,i\/E^is 
and missing momentum Pmiss = I ZliPil were calculated. To charged particle tracks the 
pion mass was assigned while the mass of clusters was assumed to be zero. 

Hadronic event candidates were required to pass the following selection criteria: 

• The total energy deposited in the electromagnetic calorimeter had to exceed 1.2 GeV 
(0.2 GeV) for ^ < 16 GeV, 2.0 GeV (0.4 GeV) for 16 < < 24 GeV and 3.0 GeV 
(0.4 GeV) for ^/s > 24 GeV in the barrel (each endcap) of the detector. 

• The number of good charged particle tracks was required to be greater than three 
reducing r^r^ and two-photon backgrounds to a negligible level. 

• For events with exactly four tracks configurations with three tracks in one hemi- 
sphere and one track in the opposite hemisphere were rejected. 

• At least three tracks had to have more than 24 hits in r — and a momentum larger 
than 500 MeV; these tracks are called long tracks. 

• The visible energy had to fulfill E^^s/ ^/s > 0.5. 

• The momentum balance had to fulfill pbai < 0.4. 

• The missing momentum had to fulfill Pmiss/ < 0.3. 

• The z-coordinate of the reconstructed event vertex had to lie within 15 cm of the 
interaction point. 
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• The polar angle of the thrust axis was required to satisfy | cos(^t)| < 0.8 in order 
that the events be well contained in the detector acceptance. 

The numbers of selected events for each y/s are shown in table 1 for the two versions of 
the data. 

3.5 Corrections to the Data 

All selected tracks and the electromagnetic calorimeter clusters remaining after correcting 
for double counting of energy as described above were used in the evaluation of the event 
shape moments. The values of the moments after all selection cuts had been applied are 
said to be at the detector level. 

In this analysis events from the process e+e~ bb were considered as background, 
since especially at low y/s the large mass of the b quarks and of the subsequently produced 
B hadrons will influence values of the moments. The QCD predictions were calculated 
for massless quarks and thus we chose to correct our data for the presence of bb events. 

For the determination of the moments we calculated the sums J^iUi'data {n = 1, . . . ,5) 
where i denotes all selected events. The expected contribution of bb background events 
^iVihh subtracted from the observed J2iyi,data- The effects of detector acceptance 
and resolution and of residual ISR were then accounted for by a multiphcative correction 
procedure. 

Two sets of J2i were calculated from Monte Carlo simulated signal events; the first, 
at the detector level, treated the Monte Carlo events identically to the data, while the 
second, at the hadron level, was computed using the true momenta of the stable particles 
in the event^, and was restricted to events where a/s', the centre-of-mass energy reduced 
due to ISR, statisfied a/s — Vs' < 0.15 GeV. The Monte Carlo ratio of the hadron level to 
the detector level for each moment value was used as a correction factor for the data. The 
corrected sums were then normalized by the expected total number of events yielding the 
final values of (y") . The expected total number of events was calculated from the number 
of selected events in the data in the same way as for the moments. 

The detector correction factors C^^*^^*°^ as determined using PYTHIA and HERWIG 
are shown in figures 1 and 2. We observe some disagreement between the detector correc- 
tions calculated using PYTHIA or HERWIG at low ^/s while at larger ^/s the correction 
factors agree well for most observables. The difference in detector corrections was evalu- 
ated as an experimental systematic uncertainty. 

4 Systematic Uncertainties 

Several sources of possible systematic uncertainties were studied. All systematic uncer- 
tainties were taken as symmetric. Contributions to the experimental uncertainties were 

estimated by repeating the analysis with varied cuts or procedures. For each systematic 
variation the value of as was determined and then compared to the result of the standard 
analysis (default value). For each variation the difference with respect to the default value 

•^AU charged and neutral particles with a lifetime larger than 3 x 10~^°s were considered stable. 
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was taken as a systematic uncertainty. In the cases of two-sided systematic variations the 
larger deviation from the default value was taken as the systematic uncertainty. 

• In the standard analysis the data version from 9/87 was used. As a variation a 
different data set from 5/88 was used. 

• In the default method the tracks and clusters were associated and the estimated 
energy from the tracks was subtracted. As a variation all reconstructed tracks and 
all electromagnetic clusters were used. 

• The thrust axis was required to satisfy |cos(^^t)| < 0.7. With this more stringent 
cut events were restricted to the barrel region of the detector, which provides better 
measurements of tracks and clusters compared to the endcaps. 

• Instead of using PYTHIA for the correction of detector effects as described in sec- 
tion 3.5, events generated with HERWIG were used. 

• The requirement on missing momentum was dropped or tightened to Pmiss/ Vs < 
0.25. 



• The requirement on the momentum balance was dropped or tightened to pbai < 0.3. 

• The requirement on the number of long tracks was tightened to Mong > 4. 

• The requirement on the visible energy was varied to E^^IVs > 0.45 and E^^/^/s > 
0.55. 

• The amount of subtracted bb background was varied by ±5% in order to cover 
uncertainties in the estimation of the background fraction in the data. 

All contributions listed above were added in quadrature and the result quoted as the 
experimental systematic uncertainty. The dominating effects were the use of the different 
data version followed by employing HERWIG to determine the detector corrections. In 
the fits of the QCD predictions to the data two further systematic uncertainties were 
evaluated: 



• The uncertainties associated with the hadronization correction (see section 5.2) 
were assessed by using HERWIG and ARIADNE instead of PYTHIA. The larger 
change in cts resulting from these alternatives was taken to define the hadronization 
systematic uncertainty. 

• The theoretical uncertainty, associated with missing higher order terms in the theo- 
retical prediction, was assessed by varying the renormalization scale factor x^. The 
predictions of an all-orders QCD calculation would be independent of x^, but a 

finite order calculation such as that used here retains some dependence on x^. The 
renormalization scale x^ was set to 0.5 and 2. The larger deviation from the default 
value was taken as theoretical systematic uncertainty. 
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5 Results 



5.1 Values of Event Shape Moments 

The first five event shape moments for the six observables after subtraction of bb back- 
ground and correction for detector effects are shown in figures 3 and 4. Superimposed 
we show the moments predicted by the PYTHf A, HERWfG and ARfADNE Monte Carlo 
models tuned by OPAL to LEP 1 data. In order to make a more clear comparison be- 
tween data and models the lower plots show the differences between data and each model 
divided by the combined statistical and experimental error for ^/s = 14 and 35 GeV. The 
three models are seen to describe the data fairly well; PYTHIA and ARIADNE are found 
to agree better with the data than HERWIG. 

5.2 Determination of cks 

Our measurement of the strong coupling constant as is based on the fits of QCD pre- 
dictions to the corrected moment values, i.e. the data shown in figures 3 and 4. The 
theoretical predictions using the 0{a\) calculation as described in section 2 provide val- 
ues at the parton level. In order to confront the theory with the hadron level data, it is 
necessary to correct for hadronization effects. The moments were calculated at hadron 
and parton level using PYTHIA and, as a cross check, with the HERWIG and ARIADNE 
models. The data points were then multiplied by the ratio C^'^'^ of the parton and hadron 
level moment values in order to correct for hadronization. 

The models use cuts on quantities like e.g. the invariant mass between partons in 
order to regulate divergencies in the predictions for the parton shower evolution. As a 
consequence in some events no parton shower is simulated and the orginal quark-antiquark 
pair enters the hadronization stage of the model directly. This leads to a bias in the 
calculation of moments at the parton level, since y = in this case for all observables 
considered here. In order to avoid this bias we excluded in the simulation at the parton 
level events without at least one radiated gluon. At the hadron and detector level all 
events were used. 

The hadronization correction factors as obtained from the three models are shown 
in figures 5 and 6. We find that the hadronization corrections reach values of down to 
about 0.5 at low ^/s. For larger ^/s the hadronization corrections decrease as expected. 
We also observe that the models don't agree well for moments based on Sw, ^23 and Mh 
at low ^/s. The differences between the models were studied as a systematic uncertainty 
in the fits. 

A x^-value for each moment (y") was calculated using the following formula: 

= - (y")theo,JV^.' (3) 

i 

where i denotes the energy points. The value was minimized with respect to q;s(-^zo) 
for each moment n separately. The running of as was implemented in the fit in two- 
loop precision using the formula shown in [26]. The scale parameter x^, as discussed in 
section 2, was set to 1. 
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The fit results are shown in figure 7. The fit to the first moment (Mh) did not converge 
and therefore no result is shown. We observe values of x^/d.o.f. of 0{1); the fitted QCD 
predictions including the running of «s ^'I'c thus consistent with our data. However, we 
find that for ((1 — T)"), (C") and (-Bt) the fitted values of as{Mzo) increase steeply with 
the order n of the moment used. This effect is not as pronounced for (-Bw), ((2/23)") 
and (Mg), n — 2 . . .5. In order to investigate the origin of this behaviour we show in 
figure 8 the ratio K = Bn/An of NLO and LO coefficients for the six observables used 
in our fits. There is a clear correlation between the steeply increasing values of as{Mzo) 
with moment n and increasing values of K with n for ((1 — T)"), (C") and (B^). The 
other observables {B^), ((1/23)") and (Mg), n — 2. . .5, have fairly constant values of K 
and correspondingly stable results for as{Mzo). We also note that (Mh) has a large and 
negative value of K which is the cause that the fit did not converge. 

It is also of interest to combine the measurements of 0:3 (Mzo) from the various fits in 
order to determine a single value. This problem has been subject of extensive study by 
the LEP QCD working group [27], and we adopt their procedure here. 

In brief the method is as follows. The measurements of as{Mzo) were combined in a 
weighted mean, to minimize the between the combined value and the measurements. If 
the measured values of 0:3 (Mzo) are denoted as,i, with covariance matrix V, the combined 
value, as (Mzo), is given by 

q;s(Mzo) = ^Wi^s.i where Wi^—^——^, (4) 

where i and j denote the individual results. The difficulty resides in making a reliable 
estimate of V in the presence of dominant and highly correlated systematic errors. Small 

uncertainties in the estimation of these correlations can cause undesirable features such as 
negative weights. For this reason only statistical correlations and experimental systematic 
errors assumed to be partially correlated between measurements were taken to contribute 
to the off-diagonal elements of the covariance matrix: Vl^ = min^al^p ^, a^^pj) . All error 
contributions (statistical, experimental, hadronization and scale uncertainty) were taken 
to contribute to the diagonal elements. The hadronization and scale uncertainties were 
computed by combining the 0:3 (Mzo) values obtained with the alternative hadronization 
models, and from setting = 0.5 or = 2.0, using the weights derived from the 
covariance matrix V. 

We considered only those results for which the NLO term in equ. (2) is less than half 
the LO term (i.e. \Kas/2n\ < 0.5 or \K\ < 25), namely (1 - T), (C), (Bt), (B^) and 
((1/23)"), n = 1, ... ,5 and (M^), n = 2, . . . , 5; i.e. results from 17 observables in total. 
The purpose of this requirement was to select observables with an apparently converging 
perturbative prediction. The K values are shown in figure 8. The statistical correlations 
between the 17 results were determined using Monte Carlo simulation at the hadron level. 

The result of the combination is 

ag(Mzo) = 0.1286 ± 0.0007(stat.) ± 0.0011(exp.) ± 0.0022(had.) ± 0.0068(theo.) , 

above but still consistent with the world average value of as(Mzo) = 0.1182±0.0027 [28]. 
It has been observed previously in comparisons of distributions of event shape observables 
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with NLO QCD predictions with — 1 that fitted values of as{Mzo) tend to be large, 
see e.g. [29]. 

Combining only the fit results from (1 — T), (C), (i^x), {Bw), (z/23) and (M^) yields 
a value of 

as{Mzo) = 0.1239 ± 0.0004(stat.) ± 0.0008(exp.) ± 0.0009(had.) ± 0.0059(theo.) . 

The slightly smaller error of as reflects the fact that the lower order moments are less 
sensitive to the multijet region of the event shape distributions. This leads to a smaller 
systematic uncertainty. 

6 Summary 

In this note we present preliminary measurements of the strong coupling constant from 
moments of event shape distributions at centre-of-mass energies between 14 and 44 GeV 
using data of the JADE experiment. The predictions of the PYTHIA, HERWIG and 
ARIADNE Monte Carlo models tuned by OPAL to LEP 1 data are found to be in rea- 
sonable agreement with the measured moments. 

From a fit of ^(as) predictions to selected event shape moments corrected for exper- 
imental and hadronization effects we have determined the strong coupling constant to 
be as{Mzo) — 0.1286 ± 0.0072 (total error). The higher moments, in particular of the 
Thrust, C-Parameter and St event shape, lead to systematically enlarged value of as. 

References 

[1] B. Naroska: Phys. Rep. 148 (1987) 67 

[2] JADE Coll., P.A. Movilla Fernandez, O. Bicbel, S. Bethke, S. Kluth, P. Pfeifen- 
schneider et al.: Eur. Phys. J. C 1 (1998) 461 

[3] JADE and OPAL Coll., P. Pfeifenschneider et al: Eur. Phys. J. C 17 (2000) 19 

[4] P.A. Movilla Fernandez: In: ICHEP 2002: Parallel Sessions, S. Bentvelsen, 
P. de Jong, J. Koch, E. Laenen (eds.), 361. North- Holland, 2003 

[5] S. Brandt, Ch. Peyrou, R. Sosnowski, A. Wroblewski: Phys. Lett. 12 (1964) 57 

[6] E. Fahri: Phys. Rev. Lett. 39 (1977) 1587 

[7] R.K. Ellis, D.A. Ross, A.E. Terrano: Nucl. Phys. B 178 (1981) 421 

[8] T. Chandramohan, L. Clavelh: Nucl. Phys. B 184 (1981) 365 

[9] L. Clavelh, D. Wyler: Phys. Lett. B 103 (1981) 383 

[10] S. Catani, G. Turnock, B.R. Webber: Phys. Lett. B 295 (1992) 269 

[11] S. Catani et al.: Phys. Lett. B 269 (1991) 432 



11 



[12] S. Catani, M.H. Seymour: Phys. Lett. B 378 (1996) 287 

[13] R.K. Ellis, W.J. Stirling, B.R. Webber: QCD and Collider Physics. Vol. 8 of Cam- 
bridge Monographs on Particle Physics, Nuclear Physics and Cosmology, Cambridge 
University Press (1996) 

[14] P.A. Movilla Fernandez: Ph.D. thesis, RWTH Aachen, 2003, PITHA 03/01 

[15] T. Sjostrand: Comput. Phys. Commun. 82 (1994) 74 

[16] G. Marchesini et al.: Comput. Phys. Commun. 67 (1992) 465 

[17] G. Corcella et al.: J. High Energy Phys. 01 (2001) 010 

[18] E. Elsen: Detector Monte Carlos, JADE Computer Note 54. 

[19] E. Elsen: Multihadronerzeugung in e'^e' Vernichtung bei PETRA-Energien und Ver- 
gleich mit Aussagen der Quantenchromodynamik, PhD thesis, Universitat Hamburg, 
1981. 

[20] C. Bowdery and J. Olsen: The JADE SUPERVISOR Program, JADE Computer 
Note 73. 

[21] L. Lonnblad: Comput. Phys. Commun. 71 (1992) 15 

[22] R. Kuhn, F. Krauss, B. Ivanyi, G. Soff: Comput. Phys. Commun. 134 (2001) 223 

[23] JADE Coll., W. Bartel et al.: Phys. Lett. B 88 (1979) 171 

[24] JADE Coll., W. Bartel et al.: Phys. Lett. B 129 (1983) 145 

[25] JADE Coll., S. Bethke et al.: Phys. Lett. B 213 (1988) 235 

[26] Particle Data Group, S. Eidelman et al: Phys. Lett. B 592 (2004) 1 

[27] The LEP Experiments (ALEPH, DELPHI, L3 and OPAL) and the LEP QCD Work- 
ing Group: paper in preparation 

[28] S. Bethke: MPP-2004-88, hep-ex/0407021 (2004) 

[29] OPAL Coll., P.D. Acton et al: Z. Phys. C 59 (1993) 1 



12 



Figures 




n 



n 




- 1 1 


1 1 


1 - 


- Pythia 






— Herwig 

- 1 1 


n 

<B^ > 

1 1 


1 - 




n 



n 



Figure 1: The figures show the detector corrections at ^/s — 14 GeV as calculated using 
PYTHIA and HERWIG. 
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Figure 2: Same as figure 1 for ^/s — 35 GeV 
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Figure 3: The figures show the first five moments of 1 — T, Mh and C at hadron level 
compared with predictions based on PYTHIA, HERWIG and ARIADNE Monte Carlo 
events. The errors shown include all statistical and experimental uncertainties. The 
lower panels show the differences between data and Monte Carlo at ^/s — 14 and 35 GeV, 
divided by the quadratic sum of the statistical and experimental error. 
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Figure 4: Same as figure 3 for B^, i^w and 7/23- 
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Figure 5: The figures show the hadronization corrections at -s/s = 14 GeV as calculated 
using PYTHIA, HERWIG and ARIADNE. 
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Figure 6: Same as fig. 5 for ^/s — 35 GeV. 
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Figure 7: Measurements of as{Mzo) using fits to moments of six event shape observables. 
The inner error bars represent statistical errors, the middle error bars include experimental 
errors and the outer error bars show the total errors. The dotted line indicates the 
weighted average described in the text; only the measurements indicated by solid symbols 
were used for this purpose. 
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Figure 8: The ratio K = Bn/An of NLO and LO coefficients for the six observables. 
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